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A complete global analysis of spatially-flat, four-dimensional cosmologies derived from the type 
IIA string and M-theory effective actions is presented. A non-trivial Ramond-Ramond sector is 
included. The governing equations are written as a dynamical system. Asymptotically, the form 
fields are dynamically negligible, but play a crucial role in determining the possible intermediate 
behaviour of the solutions (i.e. the nature of the equilibrium points). The only past-attracting 
solution (source in the system) may be interpreted in the eleven-dimensional setting in terms of flat 
' space. This source is unstable to the introduction of spatial curvature. 
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■ I. INTRODUCTION 

O ' 



There are five anomaly-free, perturbative superstring theories [Q. It is now widely believed that these theories 
represent special points in the moduli space of a more fundamental, non-perturbative theory known as M-theory p| . 
(For a review see, e.g., Ref. ||). Moreover, another point of this moduli space corresponds to eleven-dimensional 
Q\ \ supergravity. This represents the low-energy limit of M-theory ||,|| . 

The original formulation of M theory was given in terms of the strong coupling limit of the type IIA superstring. 
In this limit, an extra compact dimension becomes apparent, with a radius, R, related to the string coupling, g s , 
by R oa gl 0- The compactification of M-theory on a circle, S 1 , then leads to the type IIA superstring. In this 
framework, the dilaton field of the ten-dimensional string theory is interpreted as a modulus field parametrizing the 
,—, • radius of the eleventh dimension. 

This change of viewpoint re-establishes the importance of eleven-dimensional supergravity in cosmology and has 
interesting consequences for the dynamics of the very early universe. An investigation into the different cosmological 
models that can arise in M-theory is therefore important and a number of solutions to the effective action have 
recently been found ]|-|| . 

The bosonic sector of eleven-dimensional supergravity consists of a graviton and an antisymmetric, three-form 
potential |||. The purpose of the present paper is to employ the theory of dynamical systems to determine the 
qualitative behaviour of a wide class of four-dimensional cosmologies derived from this supergravity theory. We 
compactify the theory to four dimensions under the assumption that the geometry of the universe is given by the 
product M4 x Y 6 x S 1 , where M4 is the four-dimensional spacetime, Y 6 represents a six-dimensional, Ricci-flat 
internal space and S 1 is a circle corresponding to the eleventh dimension. We assume that the only non-trivial 
components of the field strength of the three-form potential are those on the M4 x S 1 subspace. 

The outline of the paper is as follows. In Section II we derive an effective, four-dimensional action by employing 
the duality relationship in four dimensions between a p-form and a (4 — p)-form. The field equations for the class 
of spatially isotropic and homogeneous Friedmann-Robertson-Walker (FRW) universes are derived in Section III 
and expressed as a compact autonomous system of ordinary differential equations. All of the equilibrium points of 
the system and their stability are determined in Section IV. A complete analysis of the fiat cosmological models is 
presented in Section V together with a discussion and interpretation of the results. The robustness of the models 
is addressed in Section VI wherein a number of generalizations (i.e. additional degrees of freedom) are included; in 
particular curvature effects are considered. We conclude with a discussion in Section VII. 



1 



II. FOUR-DIMENSIONAL EFFECTIVE ACTION 



The bosonic sector of the effective supergravity action for the low-energy limit of M-theory is given in component 
form byQ 



S = Id 11 
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(i) 



where R is the Ricci curvature scalar of the eleven-dimensional manifold with metric g^ v , g = det^^ and F^iA^A Aa = 
49r^ i J ^-A 2 A 3 A i ] * s ^ ne four-form field strength of the antisymmetric three-form potential A^^^. The topological 
Chern-Simons term arises as a necessary consequence of supersymmetry |J . 

In deriving a four-dimensional effective action from Eq. (pi) , we first consider the Kaluza-Klein compactification on 
a circle, S 1 . This results in the effective action for the massless type IIA superstring P,|LT|. The three-form potential 
A abc reduces t° a three-form potential Aabc an d a two- form potential, Bab = Aabt- If we ignore the one-form 
potential that arises from the dimensional reduction of the metric, the ten-dimensional action is given by fill 



S = J d w x^/\gZ\ 



( R s + (V<J> 10 ) 2 - ^H ABC H ABC 



x TP ipABCD 

-^FabcdF — 
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(2) 



where Habc = 39[aBbc] an d Fabcd — ^d[AA BC D] are the field strengths of the potentials Bsc and Abcd, 
respectively, the ten-dimensional dilaton field, $10, is related to the radius of the eleventh dimension, e 7 M: 



7 = 3*10 

and we have performed a conformal transformation to the string frame: 

g A % = n 2 g AB , il 2 = e\ 



(3) 



(4) 



The first line in Eq. (g) contains the massless excitations arising in the Neveu-Schwarz/Neveu-Schwarz (NS-NS) 
sector of the type IIA superstring and the second line is the Ramond-Ramond (RR) sector of this theory (!]]. In 
general, the NS-NS fields couple directly to the dilaton field in the string frame, but the RR fields do not. 

We now consider the compactification of theory @ to four dimensions. The simplest compactification that can be 
considered is on an isotropic six-torus, where the only dynamical degree of freedom is the modulus field parametrizing 
the volume of the internal space. We therefore assume that the string-frame metric (Q) has the form 

ds 2 s = g$dx^dx u + e 2ff S ij dx i dx j , (5) 

where Sij (i, j = 1, . . . , 6) is the six-dimensional Kronecker delta and f3 represents the modulus field. 

Moreover, we compactify the form-fields in Eq. (Q) by assuming that the only non-trivial components that remain 
after the compactification are those associated with the external spacetime M4. This implies, in particular, that the 
Chern-Simons term is unimportant, since it is proportional to F A F. The effective four-dimensional action is then 
given in the string frame by 



In this paper, the spacetime metric has signature (—,+,■■■,+) and variables in eleven dimensions are represented with a 
circumflex accent. Upper case, Latin indices with circumflex accents take values in the range A — (0, 1, ... , 10), upper case, 
Latin indices without a circumflex accent vary from A = (0,1, ... , 9), lower case Greek indices span fi — (0, 1, 2, 3) and lower case 
Latin indices represent spatial dimensions. A totally antisymmetric p-form is defined by A p = (l/p\)AA 1 ...A p dx Al A ■ ■ ■ Adx Ap 
and the corresponding field strength is given by F p +i — dA p = [l/(p + 1)!]-Fai...a +1 dx Al A ■ ■ ■ A dx Ap+1 . The coordinate of 
the eleventh dimension is denoted by T. The eleven-dimensional Planck mass is the only dimensional parameter in this theory 
fnOl and units are chosen such that 16nG = 1. 
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R+ (V$) 2 - 6(V/3) 2 
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(0) 



where 



6/3 



(7) 



is the four-dimensional dilaton field. 

The field equations and Bianchi identities for the form fields are 

V M (e-*i?^ A ) = 
d^HuXn] = 



and 



V M ^F^ Xk ) = 
9\pF u x Kp ] = 0, 



(8) 
(9) 



(10) 
(11) 



resp ectively. In four dimensions, a p-form is dual to a (4 — p)-form and Eqs. (JSj) and ( |10| ) are solved by the ansatze 



(12) 
(13) 



where e^ AK is the covariantly constant four-form, a is a scalar variable and Q is an arbitrary constant. Although 
Eqs. (0) and (|l^) solve the field equations (|J) and (|l0|), the Bianchi identities (||) and (|TT|) m ust also be satisfied. 
Eq. (|13|) is trivially satisfied, since we are working in four dimensions and substituting Eq. (|12|) into Eq. (||) implies 
that 



V M (e* W) = 0. 



(14) 



Eq. (0) may be interpreted as the field equation for the pseudo-scalar axion field, a f0| . Moreover, substituting 
Eqs. (|12|) and (|TJ|) into the remaining field equations for the graviton, dilaton and modulus fields implies that they 
may be derived from a dual effective action 
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R + (V$) 2 - 6 (V/3) 2 - ie 2 * (Vct^ 2 



(15) 



In the following Section, we derive the cosmological field equations from the effective action (|TJ|) . 



III. COSMOLOGICAL FIELD EQUATIONS 



We denote the FRW metric on M4 by the line element ds 2 



e 2a <Kl\, where e Q represents the scale 



factor of the universe, n is the lapse function and dfl 2 is the three-metric on the surfaces of isotropy, with positive 
(k = +1), negative (k = — 1) or zero (fc = 0) curvature, respectively. Substituting this ansatz into the effective action 
(0) , integrating over the spatial variables and normalizing the comoving volume to unity, yields the reduced action 



S= dt 



3d 2 ^ 2 



6kn 2 e 2a 



•6/3' 



where a dot denotes differentiation with respect to t and 

ip = $ - 3a 

defines the 'shifted' dilaton field Q. 
The corresponding field equations are 



1 

2n' 



e* +6a & 2 - 



)2 e -6/3+3a 



(16) 



(17) 
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a = a<p + ^p - 2ke~ 2a - iQV 6/3+v+3Q , (18a) 

(p = ^ \3a 2 + p 2 + 6ke- 2a + 6/3 2 - ip) , (18b) 

p = — 6ap , (18c) 

(3 = f3p + -Q 2 e - 6 ^+ 3a , (18d) 

= 3a 2 - Lp 2 + 6/3 2 + \p- 6kc- 2a + ^Q^ e ~ef3+ v +3a _ ( 18e ) 



where we specify n = 1 and 



p = e 2lfi+6a & 2 (19) 



parametrizes the kinetic energy of the pseudo-scalar axion field. 

Kaloper, Kogan and Olive have considered the equivalent compactification of the M-theory effective action (|l|) 
directly in terms of eleven-dimensional variables when M4 corresponds to the spatially flat FRW spacetime 0. In 
this case, the only non-trivial components of the four-form field strength that can exist entirely on the subspace 
A/4 x S* 1 are Fo mnp and Fxmnp, where m = (1, 2, 3), etc. The former represent the non-trivial components of the RR 
four-form field strength in Eq. ( |l5| ) and the latter are equivalent to those of the NS-NS three-form field strength. 
The scale factors of the universe in the string and M-theory interpretations are related by Eqs. (^|) and (|J). These 
relationships provide the recipe that allows the type IIA string cosmologies to be reinterpr eted in term s of eleven- 
dimensional, M-theory models. It can be verified by direct comparison that for k = 0, Eqs. (|l8a| )— (|l8e| ) are formally 
equivalent to the field equations derived in Ref. W. The advantage of employing the string- frame variables in this 
work is that the first derivative of the shifted dilaton field ([l7]) is a dominant variable and this greatly simplifies the 
analysis of the global dynamics. 



To proceed, we define a new time variable, 77: 



dq 
~dt 



3 (-6/3+ v +3q)/2 



The system of equations ( |lg| ) then becomes: 



a" = i a y + 3a'/?' 


-2 {a) 




6(/3') 2 






<p" = 


3(a') 2 + 


6 OS')' + 




\ (p'f + 3/3V - 








(3" = 




f3(/3') 2 -^'/3' 




I rtr -(v+3a-6/3) 

2 P 




3(af- 


6(/3') 2 


+ 6fce- (5Q+¥> ~ 6/3) 


2^ 



(20) 

(21) 
(22) 
(23) 
(24) 



where a prime denotes differentiation with respect to ry and the Hamiltonian constraint (18e) has been employed to 
eliminate the axion field, p. 

Since Q 2 is semi-positive definite, it follows from Eq. (24) that p' is a dominant variable in the spatially flat and 
negatively curved models (fc < 0). In addition, it follows from Eqs. (22) and (24) that ip' is a monotone function. This 
is important because it implies the global result that <p' is either monotonically increasing or decreasing throughout 
the evolution of the models. (The variable tp' plays an analogous role to that of the expansion parameter in the 
spatially homogeneous perfect fluid models of general relativity flltl ). 

We therefore introduce the new dimensionless time variable, r, according to 

dr , , . 

*T ^ ■ 

where we assume here that ip' > 0. (The case ip' < is discussed in Section V). We also define the following set of 
dimensionless variables: 



V3a' _ V6/3' _ Q 2 _ 6 fce-( 5Q +^ 6 « _ 



p< ' a " p> ' " 2(^') 2 ' " (^) 2 ' " 2(c^ 2 



n = ^2 ■ (26) 
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This leads to a decoupling of the equation for ip' , which can be written as 

a/3 , 2 2 , a/6 



dtp' 
dr 



\ + \z-—x + x'+y* + —y\ V >' 



(27) 



The remaining equations can then be written in the following dimensionless form: 



dx 
dr 
dy 
dr 

dz 

Ih 
dv 

d7 



1 / \ 2 

(l — x 2 — y 2 — z) (x + a/3) + — z yx — a/3J — ~^J^ U ; 

(l - x 2 - y 2 - z) y + (y + a/6 

- 1 - a/% + V3x + 2 (l - x 2 - y 2 - z) 



— ii 
3 



The variable f2 is given by 



2v / 3a; + 3(z + 2a; 2 + 2y 2 ) 



Q = 1 — x 2 — y — z — u, 



and satisfies the auxiliary equation 



dn 
17 



2y 2 + z + 2x (x + a/3^ 



(28a) 
(28b) 
(28c) 
(28d) 

(29) 
(30) 



IV. STRUCTURE OF STATE SPACE AND LOCAL ANALYSIS 



In this Section we present all of the equilibrium points that arise in the system (|28|). We are primarily interested 
in the spatially flat models. However, we also consider the stability of these models to perturbations in the spatial 
curvature. The local stability analysis we perform is valid for both positive and negative spatial curvature, although 
we explicitly consider the k < models since in this case the condition p > implies that all of the dimensionless 
variables are bounded. The physical state space is defined by 

< {x 2 ,y 2 ,z,u} < 1 (31) 

and a global analysis can therefore be undertaken. 

We include the boundary of the state space in our analysis because the dynamics in the invariant boundary 
submanifolds is useful in determining the global properties of the orbits in the physical phase space. The boundary 
of the state space consists of a number of invariant submanifolds of the system. They are: (i) models where the axion 
field is trivial (£1 — 0), (ii) the spatially flat models (u — 0), and (iii) models where z ~ 0, corresponding to the 
case where the four-form field strength, F^ v \ Kl is dynamically unimportant. The system of equations also admits an 
invariant submanifold, K, that is not part of the boundary of state space: 

K : x + V2y + \/3 = , u = . (32) 

The equilibrium points are: 
Equilibrium set: The Line L + 

a; 2 + y 2 = l,z = 0,u = ; O = . 



Ai = , A 2 = -2(\/3x + 1) , A 3 = -1 + \/3x - \/6?y , A 4 = -2 (-^=x + 1 

where denote the eigenvalues. The zero eigenvalue indicates that this is indeed an equilibrium set, corresponding 
to a circle of unit radius in the (x,y) plane. We refrain from presenting the eigenvectors, but note that it is the 
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eigendirection associated with A3 that points in a direction outside the submanifold z = 0, while the eigenvector of 
A4 extends into the u direction. The stability of these equilibria is discussed in the following Sectionn. 



Equilibrium point: The source R 



5\/3 V6 28 252 n „ 

, v = , z — , u = : si = 

19 19 361 361 



Al - 19 ' Aa ~ 19 



A3, 4 — 



7±n/119 



19 



Equilibrium point: The saddle M 



V3 ' 



, u 




The saddle point M corresponds to the Milne form of flat space. This may be mapped onto the future light cone of 
the origin of Minkowski spacetime and in this sense may be interpreted as the string perturbative vacuum represented 
in terms of non-standard coordinates. 



V. DYNAMICS OF THE SPATIALLY FLAT COSMOLOGIES 



A. Global Analysis 



In this Section we consider the global dynamics of the spatially flat cosmologies (k = 0, U = 0). For these models, 
the state space is three-dimensional and the orbits can therefore be represented pictorially. 

The only equilibrium points in the spatially flat models lie on the line L + and the eigenvalues are given in Section 
IV. From these eigenvalues, it can be seen that this line is a sink for x > — 1/\/3 and y/2y > x — l/y/3. The lines 
x = — 1/\/3 and \/2y = x — 1/V3 intersect on L + at the point P : (x,y) = (— 1/\/3, —y/2/3), at which all three 
eigenvalues are zero. Hence, P is a non-hyperbolic equilibrium point. All other points on L + are saddles. 

It can be shown that the point P is a source in the three-dimensional phase space. It follows from Eqs. (28a) and 
( H| ) that 



— (x + V2y + VTj = I r 



V2y + V3 ) ( 1 - x 2 - y 2 - -z 



(33) 



for u = 0. This implies that x + \/2y + v3 is a monotonically increasing function in the physical phase space. The term 
(1 — x 2 — y 2 — iz) is positive-dehnite in the interior region and can only be zero on the boundary, where x 2 + y 2 = 1 

and z = 0. The term x + \/2y + V3 is positive-definite in the physical state space and can only vanish in the extended 
phase space at the point P. Indeed, the line x + \[2y = — v3 is tangent to the unit circle x 2 + y 2 — 1 and z = 
and actually touches it at the point P. We may conclude, therefore, that the non-hyperbolic equilibrium point P is 
indeed a source for the three-dimensional system. We have verified this by analysing the equilibrium point P using 
spherical polar coordinates and by numerical calculations. 

The dynamics on the boundary of the state space is also important when interpreting the behaviour of the orbits. 
The boundary consists of the two invariant submanifolds £1 = and z = 0. The f2 = (trivial axion field) submanifold 
can be solved analytically in terms of the variables of the state space and the solution is given by 



2 For hyperbolic equilibrium points the stability is determined by the signs of the real parts of the associated eigenvalues; in 
the case of a source (past attractor) all are positive and in the case of a sink (future attractor) all are negative. Otherwise, the 
point is a saddle. If the real part of any of the eigenvalues of an isolated equilibrium point is zero, it is non-hyperbolic and the 
stability cannot be determined directly from the eigenvalues. 
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y = -V6+ iyo + ^ ){x -^ ) , (34) 
x - V3 

where (xo,yo) represents the initial point of the orbit. Thus, orbits follow straight line paths in the (x,y) plane. 
Moreover, since by definition x < this variable is a monotonically decreasing function on this submanifold and 
y is a monotonically increasing function. The line L + is a source for x > — l/v3 and V2y < x — l\/3 and a sink 
otherwise. 

The boundary z = describes models where the four-form field strength is dynamically negligible. This submanifold 
can also be solved exactly and the orbits follow the straight line paths: 

y=^#, (35) 
x Q + Vo 

where (xq, yo) again represents the initial point of the orbit. In this case, the function x is monotonically increasing 
on this submanifold. The line L + is a source for x < — 1/V3 and a sink otherwise. 

The time-reversed dynamics of the ip' > models wc have considered thus far is equivalent to the dynamics of the 
case where ip' < 0. This follows by redefining the time variable, r: 

^ = V , (36) 
an 

so that r\ and r are both increasing or both decreasing together. If we define the other state variables as in Eq. 
(29), the variables x and y for ip 1 < are now the reflections of the variables x and y for ip' > 0, i.e., x — > —x and 
y — > — y. With the new time variable (^6|), the evolution equations (|2|) will have an 'overall' change in sign, i.e., 
dx/dr — > —dx/dr, etc. Thus, the equilibrium points arc identical in both cases, but the eigenvalues have opposite 
signs. Consequently, the dynamics of the ip' < models is the time reversal of the ip' > models, where contracting 
models for ip' > are expanding models for ip' < 0, and vice versa. 



B. Physical Interpretation 

The phase space for the spatially flat models is depicted in Figs. 1-4. Figs. 1 and 2 correspond to the invariant 
submanifolds z = and O = 0, respectively. Figs. 3 and 4 represent views of two typical orbits in the full three- 
dimensional phase space. 




FIG. 1. Phase portrait of the invariant submanifold z = 0, corresponding to the case where the RR four-form field strength is 
trivial and the NS-NS three-form field strength is dynamically important. The line L + represents a line of equilibrium points. 
Large black dots denote repellers (sources) while grey-filled dots denote attractors (sinks). The point P represents a source in 
both the two-dimensional and three-dimensional phase spaces. 
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FIG. 3. Phase diagram of the spatially flat cosmologies when both NS-NS and RR form fields are non-trivial. Note that L + 
represents a line of equilibrium points. The trajectories in Figures 1 and 2 are depicted in grey in this figure along z — and 
z = 1 — x 2 — y 2 , respectively. Small black dotes represent saddle points. See also caption to Figure 1. 
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FIG. 4. An alternative view of a different trajectory in the spatially flat phase space. See also captions to Figures 1 and 3. 

The equilibrium set L + represents solutions where the form-helds are trivial and only the dilaton and moduli fields 
are dynamically important. These are known as the 'dilaton-moduli-vacuum' solutions and have an analytical form 
given by 



p a — <=<*» \ + \±h. 

e 



e a = e a *\t\ z 

« _ $» ij.i±3ft» — 1 



/9 = e /».| t |± e V(i=3«57« ) (37) 



where {a*, <J>*,/3*, h*} are constants, and e = ±1. Note that the "— " solutions in Eq. (37), which are represented by 
the line L + , correspond to t < and, in the time-reverse case (<fi < 0), the "+" solutions of Eq. ( j37|) correspond to 
t > 0. 

Let us first consider the dynamics in the invariant submanifold 2 = 0, where the NS-NS axion field is non-trivial and 
the RR four-form field strength vanishes (Fig. 1). These trajectories represent the 'dilaton-moduli-axion' solutions 
discussed in Ref. (l7). The trajectory along y — corresponds to the solution where the internal dimensions are 
static. In this case, the universe is initially contracting (x < 0), but ultimately bounces into an expansionary phase 
(x > 0). The bounce is induced by the two-form potential. It follows from Eq. (|l^) that the field strength of this 
antisymmetric tensor field is directly proportional to the volume form of the three-space. This implies that the axion 
field may be interpreted as a membrane that is wrapped around the spatial hypersurfaces Jig] . This membrane resists 
the initial collapse of the universe and results in a bouncing cosmology. Many solutions exhibit such a bounce, but 
others collapse to zero volume. These arise when the initial kinetic energy of the modulus field (internal space) is 
sufficiently high that it can always dominate the kinetic energy of the axion field. 

In the other invariant submanifold (f2 = 0), the NS-NS two-form potential is trivial, and the RR three-form 
potential is dynamical. The cosmological constant term Q 2 in the effective action ( |l5|) may be interpreted as a 0-form 
field strength. In a certain sense, this degree of freedom plays a role analogous to that of a domain wallQ |Q- However, 
in contrast to the membrane associated with the axion field, this 'domain wall' resists the expansion of the universe. 
Thus, the majority of solutions that are initially expanding ultimately recollapse, as shown in Fig. 2. There are some 



3 In general, a solitonic p-brane is supported by the magnetic charge of a (D — p — 2)-form field strength in D spacetime 
dimensions. 







solutions where the internal space is initially evolving sufficiently rapidly that the modulus field dominates the form 
field and the expansion can proceed indefinitely. Solutions that are initially collapsing do not undergo a bounce. 

In both invariant submanifolds, the point P corresponds to an endpoint on the line of sources. In Fig. 1, the 
reflection of this point in the line y — represents the opposite end of the line of sources. This point corresponds to 
a dual solution, where the radius of the internal space is inverted. Thus, the endpoints of the line of sources in the 
invariant submanifold z = are related by a scale factor duality. 

We may now consider the dynamics in the full three-dimensional phase space, where both the NS-NS two-form 
potential and RR three-form potential are dynamically significant. Although they are asymptotically negligible, the 
interplay between these fields has important consequences. The key point is that the RR field causes the universe to 
collapse, but the NS-NS field has the opposite effect. These two fields therefore compete against one another, as can 
be seen in Figs. 3 and 4. 

The point P is the only source in the system when both form fields are present. Furthermore, it follows from 
the definitions ( p6| ) that it represents the collapsing, isotropic, ten-dimensional cosmology, where a — (3. The four- 
dimensional dilaton field, is trivial in this case. As the collapse proceeds, a typical orbit moves upwards in a 
cyclical fashion until a critical point is reached, where one of the form fields is able to dominate the dynamics. The 
orbit then shadows the corresponding trajectory in the invariant submanifold z — or CI = 0. In Fig. 3, the axion 
field dominates and causes the universe to bounce. By this time, however, the kinetic energy of the modulus field has 
become significant and the solution ultimately asymptotes to a dilaton-moduli-vacuum solution on L + . 

All sinks in this phase space correspond to solutions where the internal dimensions are expanding (y > 0). There is 
a particular point where the spatial dimensions spanning the spacetime M4 become static in the late-time limit. In 
general, however, solutions either collapse to zero volume in a finite time or superinflate (a > 0) towards a curvature 
singularity. In this sense, they correspond to pre-big bang cosmologies, since the comoving Hubble radius decreases 
p5| , p0]| . However, since the internal space is expanding, it is not clear to what extent this behaviour represents a 
realistic, four-dimensional inflationary solution. 

As discussed above, the time-reversed dynamics of the above class of models is deduced by interchanging the sources 
and sinks and reinterpreting expanding solutions in terms of contracting ones, and vice-versa. Thus, the late-time 
attractor for the time-reversed system is the expanding, isotropic, ten-dimensional cosmology located at point P. 

It is of interest to reinterpret the equilibrium points of the phase space in terms of eleven-dimensional solutions. 
Since the eleven-dimensional three-form potential is trivial on L + , these points represent 'Kasner' solutions to eleven- 
dimensional, vacuum Einstein gravity. Thus, the line L + is analogous to the Kasner ring that arises in the vacuum 
Bianchi I models of four-dimensional general relativity |ll| . 

For the compactification we have considered, the scale factors in the eleven-dimensional frame are {e^e^eT}, 
where, from Eq. (|]), a = a — 7/2 and (3 = (3 — 7/2. The 'Kasner' solutions are then given by the power laws 
a = ao/t, (3 = flo/t and 7 = 70/f, where t = J dt exp(— 7/2) and the constants of integration satisfy the constraints 
3a + 6/? + 70 = 1 and ia 2 + 6/3§ + 7 g = 1. 

These redefinitions imply that the source P corresponds to the 'Kasner' solution (aoj Ab7o) = (0,0, 1). This point 
represents the Taub form of Minkowski spacetime [16[. The relevance of this solution to the problems associated 
with the pre-big bang curvature singularity have recently been discussed J^,^|, and it is interesting from an eleven- 
dimensional point of view that such a simple solution is uniquely selected by the dynamics. The endpoints of the line 
of sinks on L + correspond to the 'Kasner' solutions (1/2,0,-1/2) and (0,2/7,-5/7), respectively, and consequently 
in both cases a subset of the scale factors are static. 

This concludes our discussion of the phase space for the spatially flat cosmologies. In the following Section, we 
consider the robustness of these models to a number of possible generalizations, including the effects of spatial 
curvature. 



VI. ROBUSTNESS OF THE MODELS 



A. Effects of Spatial Curvature 



Although the compactness of the phase space depends on the fact that k < 0, one can assume arbitrary signs for 
k in order to determine the local stability of the equilibrium points in the three-dimensional set u — with respect 
to curvature perturbations. The eigenvalue associated with u for the equilibrium points L + is always negative. This 
means that the sinks on L + (i.e., points on L + for x > — and \J~2y > x — 1/V3) remain sinks in the four- 

dimensional phase-space. In addition, this implies that the point P is now only a saddle; that is, the stability of P is 
unstable to the introduction of both positive and negative spatial curvature. 
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Since a portion of L + acts as sinks in the four-dimensional phase space, there exists the global result that the 
corresponding dilaton-moduli-vacuum solutions (|37| ) (for 3d > —dp and 6/3 > 3d — ip) will be attracting solutions 
for the spatially curved models. We may deduce further global results by restricting our attention to the negatively 
curved models (k < 0), in which case the four-dimensional phase space is compact. As discussed above, the point 
P is a only a saddle point in this extended phase space. Moreover, it follows from the analysis of Section IV that 
the only attracting equilibrium point is the point R. (There is an additional saddle M which will affect the possible 
intermediate dynamics). This source corresponds to a negatively curved model with a trivial axion field; indeed, it is 
a power-law, self-similar collapsing solution with non-negligible modulus and RR form-field. 

We have been unable to find a monotone function on the extended four-dimensional phase space, but it is plausible 
that all negatively-curved models evolve from the solution corresponding to the global source R towards the dilaton- 
moduli-vacuum solutions (on the attracting portion of L + ). Clearly the curvature is dynamically important at early 
times. 

In the time-reverse case, the solutions asymptote from the non-inflationary dilaton-moduli-vacuum solutions in the 
past and evolve to the future towards a curvature dominated model; it is plausible that they evolve towards a model 
which is the time-reversal of the one represented by R. Therefore, curvature can also be dynamically important at 
late times. 



B. Effects of Generalized Couplings 



We now consider a generalization of the effective action (15) given by 



S = J dW^sje"* 



- ■ ■ "T- r,2 1 .2* (TT7„\2 2A 



R+ (V$) z - 6(V/3) 



(V<7) 



2 e c/3 



(38) 



where A represents a cosmological constant term and c is an arbitrary constant. The former term may arise through 
non-perturbative corrections to the string effective action. The motivation for considering an arbitrary coupling of 
the modulus field to the four-form field strength is that the generality of the dynamics discussed in Section IV (in 
which c = —6) can be investigated. Eq. (|38|) reduces to the action studied in Ref. |2l| when c = 0. 
By invoking the same assumptions as in section III, the action (pq) reduces to 



S= dt 



3d — (p + 6ke' 



-la 



6/3 2 - 2A 



}2 e c/3+3a 



(39) 



and the corresponding field equations can again be derived from this action. In analogy with Eqs. (23) and (28), we 
introduce a new time variable, r, defined by 



l c -i(c/3 +¥; +3«)ii 

dt 



and the new reduced variable 



2Ae~( 3a+¥,+c/3 ) 



(40) 



(41) 



From these definitions and the reduced variables defined earlier, we obtain a five-dimensional system of ordinary 
differential equations for the reduced (dimcnsionless) variables after eliminating the variable Q that is now defined by 
£1 = 1 — x 2 — y 2 — z — u — v. Since p > 0, all of the dimcnsionless variables are bounded for the models with k < 
and A > 0, where the physical state space is defined by < {x 2 , y 2 , z, u, u} < 1, and a global analysis is therefore 
possible in this case. Including the boundaries f2 = 0, z — 0, u = and v — leads to a compact state space. 

We can analyse these models and obtain qualitative information about the dynamics in an analogous way to 
that done in section III |22|. The equilibrium set L + : 
eigenvalues associated wit 
1/V3,y = - 



1,2 



still exists, and since the 



,h[u and) v are (both) negative, part of L + will act as sinks and the non-hyperbolic point 
x = — 1/V3, y = — v2/3, z = u = v = is clearly a saddle. These are local results and are valid in all cases. 
There also exists an attracting equilibrium point W: 



V3(12 



J) D 2 



y 



8cV6D 2 , 



-192CD 4 



v = 



(42) 



with eigenvalues 



-2CD 2 , -(C± v/C(17c 2 + 924))L> 2 , -4(12 + c 2 )D 2 , -6(12 + c 2 )D 2 
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where C = (c + 6)(c — 6) and D 2 = 108 + c 2 . Note that this point corresponds to an exact self-similar collapsing 
cosmological solution with non-trivial modulus and four-form fields. The value of £1 as a function of c is given by 

-2C(60 + c 2 ) 

(108 + c 2 ) 2 ' [ 6) 

which implies that —6 < c < 6 in order for O > 0. When c = ±6, the point W is a part of the equilibrium set L + ; 
in fact, it becomes just the non-hyperbolic point P discussed previously. Note that it is also a part of the invariant 
submanifold cx — 6\/2y + \/3c = 0; u = v = 0, which generalizes the invariant submanifold K defined in Eq. (|32]). 

There are two non-flat (u ^ 0) vacuum equilibrium points with a vanishing cosmological constant v = 0, one of 
which is a source and the other a saddle. There is also a non-hyperbolic vacuum equilibrium point with a non- 
vanishing cosmological constant (v = 1) with x = y = z = u = 0, which appears to be a source. In addition, we 
can find monotone functions in the boundary submanifolds; indeed, the boundary O = u = v = and the boundary 
submanifold z = u = v = can be solved exactly in terms of the variables of the state space. Exact solutions of the 
equations of motion for particular values of c can also be found. 

However, the primary motivation for these comments is to emphasize two important points regarding the very 
interesting dynamics of the M-theory cosmologies studied earlier. First, we note that the conclusions obtained for 
the spatially curved models are robust when additional physical fields (e.g., a A term) are included. Second, and 
perhaps more importantly, we see that the value c = —6 is a bifurcation value in the analysis of general models with 
arbitrary coupling, c. In this context, therefore, the M-theory cosmological models we have studied exhibit rather 
unique dynamics. 



VII. DISCUSSION 



In this paper we have presented a complete dynamical analysis of spatially flat, four-dimensional cosmological 
models derived from the M-theory and type IIA string effective actions. We have shown that models generically 
spiral away from a source P, undergoing bounces due to the interplay between the NS-NS two-form potential and the 
RR three-form potential. Eventually, they evolve towards dilaton-moduli-vacuum solutions with trivial form fields 
(corresponding to the sinks on L + ). We note the important dynamical result that ip' is monotonia 

Thus, the form fields that arise as massless excitations in the type IIA superstring spectrum, or equivalently from 
the three-form potential of eleven-dimensional supergravity, may have important consequences in determining initial 
and final conditions in string and M-theory cosmologies, even though they are dynamically negligible in the early- 
and late-time limits. In particular, the point P is the only source in the system. It can be interpreted in the string 
context as the isotropic, ten-dimensional solution. Alternatively, it represents the Taub form of flat space when viewed 
in terms of eleven-dimensional variables. 

When the effects of spatial curvature are included, we obtained the local result that the point P becomes a saddle. 
On the other hand, the dilaton-moduli-vacuum solutions with trivial form fields are generic attracting solutions. In 
the analysis of the negatively-curved models, we found that the early time attractor (the source R) has non-zero 
curvature, implying that spatial curvature is dynamically important at early times in these examples. 

This work can be generalized in a number of ways. We considered a specific compactification from eleven to 
four dimensions, where the topology of the internal dimensions was assumed to be a product space consisting of a 
circle and an isotropic six-torus. We emphasize, however, that the analysis also applies to compactifications on a 
Calabi-Yau three-fold since the gauge fields arising from the higher-dimensional metric have been ignored pj. The 
qualitative analysis may be readily extended to compactifications on a general, rectilinear torus S 1 x . . . x S . After 
suitable redefinitions of the additional moduli fields that subsequently arise, the dimensionally reduced action can 
be expressed precisely in the form of Eq. (EH), with the inclusion of a set of massless scalar fields in the NS-NS 
sector. In particular, the compactification on T x T x S , where T n represents the isotropic n-torus, is relevant 
to compactifications involving the four-dimensional space K3 p3| . This space has played an important role in 
establishing various string dualities ||. It is the simplest four-dimensional, Ricci-flat manifold after the torus pi ] 
and may be approximated by the orbifold K3 w T 4 /Z2 pq ]. 

Moreover, the effects of spatial anisotropy in the spacetime M4 can also be considered by introducing two, uncoupled 
moduli fields into the NS-NS sector of the reduced action (|l6l). In this context, such fields parametrize the shear in 
the cosmologies. When these fields are non-trivial, the models represent the class of isotropic curvature cosmologies 
and correspond to Bianchi type I, V and IX universes [^l],^6) . It would be interesting to consider these generalizations 
further. 
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